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We present Spitzer-IRAC, NOAO 2.1 meter-Flamingos, Keck-NIRC, and FCRAO-SEQUOIA obser- 
vations of the massive star forming complex S254-S258, covering an area of 25 x 20 arc-minutes. Using 
a combination of the IRAC and NIR data, we identify and classify the young stellar objects (YSO) 
in the complex. We detect 510 sources with near or mid IR-excess, and we classify 87 Class I, and 
165 Class II sources. The YSO are found in clusters surrounded by isolated YSO in a low-density 
distributed population. The ratio of clustered to total YSO is 0.8. We identify six new clusters in 
the complex. One of them, G192. 63-00, is located around the ionizing star of the Hn region S255. 
We hypothesize that the ionizing star of S255 was formed in this cluster. We also detect a southern 
component of the cluster in Hn region S256. The cluster G192.54-0.15, located inside Hn region S254 
has a Vlsr of 17 km s" 1 with respect to the main cloud, and we conclude that it is located in the 
background of the complex. The structure of the molecular cloud is examined using 12 CO and 13 CO, 
as well as a near-IR extinction map. The main body of the molecular cloud has Vlsr between 5 and 
9 km s~ x . The arc-shaped structure of the molecular cloud, following the border of the Hn regions, 
and the high column density in the border of the Hn regions support the idea that the material has 
been swept up by the expansion of the Hn regions. 

Subject headings: Hn regions — stars: formation — stars: pre-main sequence — stars: early-type — 
infrared: stars 



1. INTRODUCTION 

The Spitzer space telescope has proved to be a useful 
tool to study the spatial distributions and evolutionary 
states of young stars in regions of star formation. The 
combination of Spitzer-IRAC and near-IR (NIR) data 
allow us to identify the young stellar objects (YSO) in 
those regions. Several studies involving IRAC+NIR data 
have appeared recently, most of them on low-mass re- 
gions and in part due to the Legacy and GTO surveys 
(Padgett et al. 2008; Chapman et al. 2007; Harvey et al. 
2007b; Porras et al. 2007; Jorgensen et al. 2006; Guter- 
muth et al. 2007; Allen et al. 2007). In this work, we 
present Spitzer-IRAC, 2.1 meter-Flamingos, and Keck- 
NIRC imaging covering the whole complex of high mass 
star formation S254-S258. We use the IRAC+NIR data 
to identify and classify the YSO in the complex. We also 
combine our data with 12 CO and 13 CO emission in the 
region, in order to understand how massive stars interact 
with the interstellar medium. In this paper, we develop 
the method of analysis and present the results for the 
S254-S258 complex. In a following paper, we will repeat 
our analysis in other massive star forming regions and 
compare the results with the studies for low mass star 
forming regions (Chavarri'a et al. 2008, in preparation). 

Following is a brief description of the S254-S258 com- 
plex. In § 2 we describe our data. Results of the analysis, 
including the identification of YSO, an analysis of their 
spatial distributions, and an investigation of the struc- 
ture of the molecular cloud, are present in § 3. In § 4 we 
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discuss individual clusters, the ages of the Hn regions, 
and the overall picture. Our conclusions are presented in 
§5. 

1.1. Sharpless 254-258 

The S254-S258 complex is a star forming region de- 
scribed for the first time by Sharpless (1959) as a group 
of five Hn regions (S254 to S258). The complex is lo- 
cated inside an area of 25x20 arc- minutes on the sky 
(see Figure 1). The ionizing stars of four of the Hn re- 
gions have been identified: a 09. 0V star in S254 (Moffat 
et al. 1979), B0V in S255 (Moffat et al. 1979), B2.5V 
in S256 (Russeil et al. 2007), and B0.5V for S257 (Reed 
2005). 

Clusters of stars have been identified in Hn regions 
S256, S258, and another inside S254 (Dutra & Bica 
2001). Until now, there is no reference to the number of 
cluster members or percentage of sources with IR-excess 
in those clusters. 

A good summary of distance estimations to the com- 
plex S254-S258 can be found in Hunter & Massey (1990) 
and Carpenter et al. (1995). In general, there is good 
agreement on a range of distances between 2 and 3 kpc. 
We calculated the average of the photometric distances 
available in the literature and adopted it as the distance 
to the complex. The photometric distances are 2.6, 2.5, 
and 2.0 kpc, from Pismis & Hasse (1976), Moffat et al. 
(1979), and Chavarria-K et al. (1987) respectively, and 
our adopted distance is 2.4 kpc. This distance is in good 
agreement with the distances estimated by Evans et al. 
(1977) (2-3 kpc, see their appendix), Pismis & Hasse 
(1982) (2.3 kpc, kinematic distance using a Schmidt ro- 
tation curve), and also Russeil et al. (2007) (2.46 kpc, 
photometric distance). 

Most of the research on the complex has been done 
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on a young massive cluster located between the Hn re- 
gions S255 and S257. This cluster, called S255-2 (also 
S255IR), is located along a dense and dusty filament 
of molecular material (Heyer et al. 1989). ft contains 
several signatures of star formation, such as infrared 
sources, OH, H2O and methanol masers, as well as 
HH-like objects (Beichman et al. 1979; Turner 1971; Lo 
& Burke 1973; Minier et al. 2005; Miralles et al. 1997). 
Approximately one arc- minute north from S255-2, there 
is a far-IR source detected by Jaffe et al. (1984), called 
S255N. This source is also associated with signposts 
of massive star formation, but since there is no NIR 
sources associated with S255N, it is believed that this 
region is in an earlier evolutionary stage than S255-2. 
Cyganowski et al. (2007) show evidence of a massive 
protocluster in S255N. They found 3 cores at 1.3 mm 
with no IR counterpart, with masses between 6-35 M . 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Mid-IR imaging 

We observed the complex S254-S258 with the Spitzer 
space telescope in December 2004 and October 2007. 
Since the IRAC field of view for the 3.6 and 5.8 ^m 
bands is shifted with respect to the field of view for 4.5 
and 8.0 /im bands, the area covered by channels 1 and 3 
is not the same as the area covered by channels 2 and 4 
(see Table 1). The IRAC pixel scale is 1.2 arc-seconds per 
pixel. We used an integration time of 10.4 seconds per 
dither in the High Dynamic Range (HDR) mode, with 3 
dithers per map position. HDR mode also acquires 0.4 
seconds integration time frames for the recovery of bright 
sources which are saturated in longer exposures. We used 
S. Carey's artifact correction scripts to remove column 
pull-down and some of the banding and muxbleed arti- 
facts (Hora et al. 2004). IRAC mosaics were constructed 
using the Basic Calibrated Data (BCD) frames (S14.0.0, 
and S15.3.0) with IRACproc (Schuster et al. 2006). The 
final image scale is 0.6 arc-seconds per pixel. 

We used IRAF DAOPHOT packages to extract sources 
and perform aperture photometry in each band. The 
photometry was done using an aperture of 1.8 arc- 
seconds for the 3.6 and 4.5 /im bands and 2.4 arc-seconds 
for the 5.8 and 8.0 ^m bands. Inner and outer sky annuli 
of 4.8 and 6 arc-seconds respectively were used in each 
IRAC band. We calculated the zero points in each band 
using Vega fluxes. The zero points used were 18.443, 
17.879, 17.234 and 16.477 for the 3.6, 4.5, 5.8 and 8.0 
/xm bands respectively. These include aperture correc- 
tions. 

2.2. Near-IR imaging 

Observations with Flamingos were performed at the 
2.1 meter telescope located at Kitt Peak National Ob- 
servatory in December 2004 and January 2006. Flamin- 
gos has a 2048 x 2048 pixel Hawaii II HgCdTe detector 
array with a plate-scale of 0.611 arc-seconds per pixel 
which gives a field of view (FOV) of 20 x 20 arc-min. We 
obtained imaging at J, H and K bands centered at 1.24 
1.65 and 2.21 ^m respectively. The seeing during the 
observations was approximately 1.0 arc-second. The ob- 
servations were done in dithering mode with 15-30 arc- 
seconds shifts for a total integration time of 1000 seconds 
per band. 



The near-IR data reduction was performed using 
IDL 4 , IRAF 5 and WCS Tools 6 routines, which includes 
linearization, bad-pixel, dark and flat-field creation and 
application, background frame creation and subtraction, 
distortion measurements and correction, mosaicking 
and astrometry correction using the positions of ap- 
proximately 2000 stars from the 2MASS catalog. Point 
source detection and aperture photometry were done 
using IRAF DAOPHOT packages. Aperture photome- 
try for saturated stars was replaced by non-saturated 
2MASS photometry and by visual inspection all the 
detections identified as nebulosity were considered as 
non-stellar and rejected. Parameters used for detection 
and aperture photometry were: threshold of 3 sigma 
for detection and 1.8, 4.8 and 6 arc-seconds in radii for 
aperture, inner annulus and inner plus outer annulus 
respectively. Calibration was performed by minimizing 
residuals to corresponding 2MASS detections. The RMS 
for residuals between the data sets used was 0.07, 0.09 
and 0.08 magnitudes for J, H and K bands respectively. 

NIRC observations of S255-2 were generously supplied 
by L. Hillenbrand. The data includes observations in H 
and K bands, that were taken using the Keck telescope 
in the year 1999. The resolution is 0.15 arc-seconds per 
pixel, with a FOV of 1 x 1 arc-minute. The observations 
were done in dithering mode with offsets between 3 and 
10 arc-seconds for a total integration time of 300 seconds 
in K-band and 400 seconds in H-band. 

We performed aperture photometry of 0.6 arc-seconds 
in radii for both bands, using a detection threshold of 
3 sigma, and 1.5 and 2.25 arc-seconds radii for inner 
annulus and inner plus outer annulus respectively. As 
for Flamingos, we calibrated NIRC photometry by min- 
imizing residuals to corresponding 2MASS detections. 

Finally, we combined IRAC, Flamingos and NIRC data 
by merging the photometry of the seven bands. In the 
area covered by both Flamingos and NIRC we used NIRC 
photometry for all the sources detected in the NIRC mo- 
saics, and Flamingos photometry for sources not detected 
in the NIRC mosaics. The maximum tolerance in posi- 
tional offsets between bands was 2.0 arc-seconds. 

2.3. Optical and NIR spectroscopy 

Optical spectra were obtained using Hectospec at the 
MMT on Mt. Hopkins. We observed the ionizing sources 
of Hn regions S254 and S255 (Reed 2005), and also a 
source in S256. Near-IR spectrum was obtained for one 
source in S258 using SpeX (Rayner et al. 2003) at the 
IRTF telescope in Mauna Kea. Hectospec data were re- 
duced using the TDC pipeline at the CfA, and classified 
using SPTclass 7 . SpeX data was reduced using Xspex- 
tool (Cushing et al. 2004), and telluric correction was 
performed using Xtellcor (Vacca et al. 2003). The sources 
in S256 and S258 were chosen because they are located 
near the center of the Hn regions, and also because they 

4 Linearization, developed by Robert A. Gutermuth. 

5 Darks, Flat-field, Bad-pixel mask and background frame cre- 
ation and application. Distortion correction and mosaicking, de- 
veloped by Luis A. Chavarria. 

6 Distortion and astrometry correction, developed by Doug 
Mink. 

7 http://www.astro.lsa.umich.edu/~hernandj/SPTclass/sptclass.html 
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are the most luminous sources inside the Hn regions at 
NIR wavelengths, and are therefore good candidates for 
being the ionizing source of the Hn regions (at optical 
wavelengths, there is one source more luminous than the 
observed source in S258, but we rejected it because it is 
located at the edge of the Hn region) . 

2.4. Millimeter 

Observations of the J=l-0 spectral lines of 12 CO and 
13 CO were conducted at the Five College Radio Astron- 
omy Observatory (FCRAO) 14 meter telescope in New 
Salem, Massachusetts. The FCRAO beam size is 45 arc- 
seconds in 12 CO and 46 arc-seconds in 13 CO. We used 
the 32 pixel SEQUOIA focal plane array (Erickson ct 
al. 1999) to image these spectral lines in the on-the-fly 
(OTF) mapping mode. The dual channel correlator al- 
lowed the two spectral lines to be observed simultane- 
ously, in 1024 channels over a bandwidth of 50 MHz 
at each frequency. The total velocity coverage exceeds 
120 km s _1 , centered on Vlsr = km s^ 1 , with chan- 
nel spacings of 0.126 km s" 1 ( 12 CO) and 0.132 km s" 1 
( 13 CO). The velocity resolution is 1.21 times the channel 
spacing. 

These data were taken as part of the Extended Outer 
Galaxy Survey (E-OGS; Brunt & Heyer, in preparation) 
that extends the coverage of the FCRAO Outer Galaxy 
Survey (OGS; Heyer et al. (1998)) to Galactic longi- 
tude / = 193, over a the latitude range -3.5< b <+5.5. 
Further details on the OTF scanning procedure can be 
found in the E-OGS data paper. During the survey ob- 
servations, pointing and focus checks were carried out 
every 3-4 hours, shortly after dawn/dusk or after a sig- 
nificant change in source coordinates. The data were 
initially converted to the T\ scale using the standard 
chopper wheel method (Kutner & Ulich 1981). We 
used the OTFTOOL software, written by M. Heyer, G. 
Narayanan, and M. Brewer, to place the spectra on a 
regular 22.5 arc-seconds grid in Galactic I, b coordinates. 
During this procedure, first order baselines were fitted 
to the spectra over velocity ranges free of spectral lines. 
These baselines were subtracted from the data, and root 
mean square noise amplitude (a) of each spectrum was 
recorded. Individual spectra contributing to a single 
Galactic coordinate were assigned a 1/er 2 weighting dur- 
ing the gridding. The gridded data were scaled to the 
main beam temperature scale by dividing by the main 
beam efficiency of 0.48. 

3. RESULTS 

The number of sources detected and the total area ob- 
served in each band is shown in Table 1. We include 
only sources with a photometric error of less than 0.2 
magnitudes in our study. The 90% completeness limit 
in each IRAC and NIR band is also shown in Table 1. 
Table 2 lists the coordinates, magnitudes and classifica- 
tion of sources. The complete table is available in the 
electronic version of this paper. Table 3 shows the result 
of our spectral classification for the ionizing stars of the 
Hn regions. 

3.1. Description of the IRAC data 

As shown in Table 1, the number of sources detected 
with IRAC decreases dramatically in channels 3 and 4 
with respect to channels 1 and 2. This difference is 



caused by the fact that channels 1 and 2 are more sen- 
sitive to the photosphcric emission from stars. Also, the 
bright extended emission in channels 3 and 4 decreases 
our sensitivity (see end of § 3.2.2). The number of sources 
detected in each of the IRAC bands per magnitude is 
shown in Figure 2. We estimate a 90% completeness at 
a magnitude of 14, 13.5, 11, and 10 for IRAC channels 1 
through 4 respectively. The completeness was estimated 
by adding and recovering artificial stars. 

The long exposure mosaic for each of the IRAC chan- 
nels is shown in Figure 3. Most of the extended emission 
seen is from polycyclic aromatic hydrocarbon (PAH) fea- 
tures in the IRAC bands (at 3.3, 6.2, 7.7, and 8.6 /nm). 
The PAH emission is detected around the five Hn re- 
gions in the complex. As shown in the 8 /im mosaic and 
in Figure 1, the PAH emission in S254 forms a shell that 
roughly follows the perimeter of the Hn region. The PAH 
emission inside the shell is very faint, with the exception 
of the young cluster located south-east from the ionizing 
star. The PAH emission in S255 and S257 has a similar 
morphology as in S254, with more PAH emission coming 
from the inner parts of the Hn regions. In S256 the PAH 
emission has a shell-like morphology with emission from 
the inner part of the Hn region. The structure of the 
PAH emission in S258 has no shell-like structure visible. 

3.2. Identifying and classifying young stars 

In this section, we describe how YSO were identified 
and classified. We classified the YSO using two meth- 
ods: a) the observed IRAC spectral energy distribution 
(SED) slope ( 

Chirac): for sources detected in all four 
IRAC bands, and b) the dereddened IRAC SED slope 
(a IRAC ), fo r sources with 7-band detections (J, H, K, and 
the four IRAC bands). We identified additional YSO by 
combining the H, K and IRAC 4.5 (im bands. Probable 
background sources were eliminated using the algorithms 
explained in § 3.2.4. 

3.2.1. Sources with IRAC 4~band detections 
The SED slope is defined as: 

dlog(XFx) m 
dlog(X) ' { ' 

and was introduced as a tool to classify YSO by Lada 
(1987). In that paper, the YSO were classified in three 
categories; Class I (star plus infalling envelope), Class II 
(star plus a disk), and Class III (post T-tauri stars). The 
classification of YSO evolved with time, adding more cat- 
egories like Class (Andre et al. 2000), or flat spectrum 
(Greene et al. 1994), as well as changing the boundaries 
between them. 

Here, we use the IRAC SED slope (a IRAC ) to classify 
the YSO using the following criteria: Class I (0 < a IRAC ), 
Class II (—2 < a IRAC < 0), and Class III or main sequence 
photospheres (a IRAC < —2). There could be some Class 
sources in our Class I category, but we are not able 
to separate them. This analysis was done for all the 
sources that were detected in all four bands. We use 
Chirac because the slope of the SED is less affected by 
extinction in the IRAC wavelengths than in the near-IR. 
To estimate a IRAC , we performed a minimum % 2 fit over 
the IRAC fluxes (see Figure 4). 

Figure 5 shows the distribution of the observed a IR AC in 
the S254-S258 complex. The distribution is similar to the 
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distribution seen in nearby star forming regions (Lada et 
al. 2006; Kumar et al. 2007; Harvey et al. 2007a), and 
also in the SED models by Robitaille et al. (2006). The 
peak corresponding to Class II sources may be explained 
by the amount of time that stars stay in a given stage of 
evolution. 

We detected a total of 462 sources in all four IRAC 
bands. Of those, using the observed a mAC , 93 were clas- 
sified as Class I and 160 as Class II. 

3.2.2. Sources with 7-band detections 

For sources that have NIR detections in J, H, and K- 
bands in addition to the four IRAC bands, we can correct 
for extinction and use the dereddened value of a mAC to 
classify them. The sources were dereddened by moving 
them to their intrinsic main sequence color in the H-K 
vs. J-H color-color diagram (Bessel & Brett 1988) us- 
ing the reddening law from Indebetouw et al. (2005) in 
JHK-bands. After the sources were dereddened, we per- 
formed a minimum x 2 fit over the dereddened fluxes in 
the IRAC wavelengths to estimate a IRAC , and then we 
classified them using the criteria explained in the previ- 
ous section. We used the extinction law from Flaherty 
et al. (2007) to deredden the sources in the IRAC wave- 
lengths. 

We detected a total of 264 sources in all 7-bands. Of 
those, using the dereddened a IRAC , 26 were classified as 
Class I and 109 as Class II. If we compare this result with 
the classification of the same sources using the observed 
Chirac: we find that 17 (6%) change their classification: 
6 Class I became Class II, and 11 Class II became Class 
III after correcting for extinction. 

Finally, by combining the two methods, that is, using 
the 4-band classification method and correcting by the 
dereddened value of a IRAC in sources with 7-band detec- 
tions, we detected in total 87 Class I, and 165 Class II out 
of 462 sources in the S254-S258 complex. Figure 6 shows 
an IRAC color-color diagram with the sources coded by 
classification. 

The PAH emission present in the complex decreases 
the sensitivity to faint sources at 5.8 and 8.0 /im in areas 
where the emission is bright. In order to quantify this 
effect, we compared the magnitudes of sources located in 
areas with bright PAH emission against sources located 
in areas with faint PAH emission. Figure 7 shows the 
distribution of YSO detected inside areas with bright 
emission (in gray), and outside those areas (white) at 
5.6 and 8.0 //m (areas with bright emission are defined 
in Fig. 3c). Also, it shows the distribution of Class I 
(in gray) and Class II (in white) sources detected at 5.6 
and 8.0 fim. From the histograms we can deduce: 1) 
the sources detected outside the areas with bright PAH 
emission extend to fainter magnitudes than inside those 
areas (this can also be seen in Figure 8), and 2) the dis- 
tribution of Class I and Class II shows that on average 
Class II sources are fainter than Class I sources. We con- 
clude that we are not detecting some of the faint YSO in 
locations with bright PAH emission, and the undetected 
YSO are preferentially Class II. The regions that are go- 
ing to be likely affected by this bias are located between 
Hn regions S255 and S257. 

3.2.3. Identifying more YSO 



We used a combination of H, K, and IRAC 4.5 
bands to detect more sources with IR-cxcess. We used 
the 4.5 /im band because it is more sensitive to IR excess 
than 3.6 fim. Since now we do not require detections in 
IRAC bands 3 and 4, we are able to detect faint sources 
not detected before due to the bright extended emission 
in those channels. 

Figure 9 shows the IRAC+NIR color-color diagram us- 
ing H, K and 4.5 /xm. In the diagram, the slope of the 
reddening vector was determined using the extinction law 
from Flaherty et al. (2007). The stars with infrared ex- 
cess are located to the right of the reddening vector and 
plotted as diamonds in the diagram. In order to clas- 
sify a source as having an IR-excess, the source must 
be located more than 1 sigma away from the reddening 
vector. 

In total, we found 395 stars with NIR excess with this 
method. Of those, 138 correspond to sources detected in 
the 4 IRAC bands; 38 Class I, 99 Class II, and 1 Class III. 

Adding all the sources with NIR excess and the Class I 
and Class II sources described in previous sections, we 
detected in total 510 sources with IR-excess in the S254- 
S258 complex. 

3.2.4. Background contamination 

Due to the sensitivity of IRAC observations and the 
relative transparency of molecular clouds at mid-infrared 
wavelengths, we expect to have field contamination in our 
sample. Results from the Spitzer c2d Survey have shown 
that galaxies have colors that are similar to the colors of 
young stars (J0rgensen et al. 2006). In order to separate 
those galaxies from the YSO, we used a color-magnitude 
analysis to isolate most (> 90%) of the background con- 
tamination. 

The data used in our background analysis were taken 
from the IRAC Shallow Survey (Eisenhardt et al. 2004). 
This survey covers 8.5 square degrees of the sky in the 
NOAO Deep Wide-Field Survey in Bootes with 30 second 
exposures, which makes it deeper than our IRAC data. 

We used two different color-magnitude diagrams in our 
analysis: [4. 5] -[8.0] vs. [4.5] for sources detected at least 
in the four IRAC bands, and K-[4.5] vs. K for the rest of 
the sources detected in K band and 4.5 /im. In order to 
compare our data with the Shallow sample, we normal- 
ized the number of galaxies by the ratio of the observed 
areas, and we added to the galaxy fluxes an extinction 
of Ak = 0.4, corresponding to the median extinction in 
the S254-S258 complex. Figure 10 shows the diagnos- 
tics used to separate galaxies from YSO for sources with 
detections in the four IRAC bands (upper) and in the 
NIR+IRAC (lower). The left part of the diagram shows 
the sources detected in the S254-S258 complex, and the 
right part shows the distribution of galaxies from the 
Shallow survey normalized by the area. We defined the 
"contaminated area" as the area with more than 1 galaxy 
per color-magnitude bin. In order to minimize the ex- 
tragalactic background, we exclude from our analysis all 
the sources detected inside the contaminated area in the 
color-magnitude diagram. 

We estimate that before correction, the number of 
galaxies in our mosaics is around 30. After applying 
these background elimination criteria, we estimate that 
there are at most 3 galaxies left in our sample. 
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3.3. Spatial Distributions of YSO 

3.3.1. Class I and Class II sources 

Figure 1 1 shows the spatial distribution of Class I and 
Class II sources, based on our classification in § 3.2. In 
general, Class II sources are more dispersed than Class I. 
There is a clear concentration of Class I sources around 
S255-2 that is extended in the north-south direction, 
forming a filament of YSO between S255 and S257. An- 
other, more dispersed group of Class I is located in S256 
and to the south of S256 and S254. We also detected 
Class I sources in S258, in the cluster near the center 
of S254, and a filamentary group in the lower part of 
the mosaic. We do not detect Class II sources in the 
core of S255-2, but this may be due to the bright PAH 
emission in the area (see end of § 3.2.2). We detected 
a high concentration of Class II sources in S256 and to 
the south between S254 and S256. There is also a halo 
of Class II sources north from S255-2, spreading to the 
north between S255 and S257. There are more Class II 
sources towards S258, a small group located south from 
S258, and another clustered group in the lower part of 
the mosaic. 

Since Class II sources are assumed to be more evolved 
than Class I, we can use the ratio of Class II to Class I 
sources as a proxy for the relative ages of the clusters. In 
this scheme, even though we are missing some Class II 
sources, the high concentration of Class I suggests that 
the youngest population is located between the Hii re- 
gions S255 and S257 (see Table 4). 

3.3.2. IR-excess sources 

Figure 12 shows the spatial distribution of all the 
sources with IR-excess (from § 3.2.3). We detect a large 
number of YSO between S255 and S257 that were not de- 
tected previously due to the extended emission at 5.8 and 
8.0 fim. The filament of YSO located between S255 and 
S257 is mixed to the south with another group of YSO 
concentrated in S256 and extended east- west to the south 
of S257 and S254. We also detected a large concentra- 
tion of YSO in S258, and other four groups: one located 
south from S258, another south-east from the first one, 
another group inside S254, and a group in the south-west 
corner of the image. 

3.3.3. IR-excess sources surface density map 

Having identified the young population of stars via 
their infrared excess, we studied the spatial distribution 
of YSO in the complex. Contours in Figure 12 corre- 
spond to the surface density map of sources with IR- 
excess. The map was created using a 3 arc-second grid. 
At each point of the grid, we calculated the surface den- 
sity of stars given by: 



where r/v is the distance to the N = 5 nearest neighbor 
(NN). The contours identify clusters S255-2, S256, S258 
and the cluster inside S254, but also show other concen- 
trations of YSO in the complex towards the south of the 
mosaic. 

Another feature seen in Figure 12 is that YSO clus- 
ters seems to be surrounded by a low-density distributed 



population of YSO. In order to determine the percent- 
age of YSO located in clusters, we need to separate them 
from the low-density component. We establish the sep- 
aration using the near neighbor cumulative distribution 
of sources with IR-excess. Figure 13 shows the NN cu- 
mulative distribution for sources with IR-excess in the 
complex S254-S258. We have chosen (somewhat arbi- 
trarily) the inflection point at 0.01 degrees (or 0.4 pc 
at 2.4 kpc of distance) as d c , the maximum separation 
between cluster members. 

Using d c , and imposing a minimum number of mem- 
bers of 5 stars per cluster and a surface density of at least 
5 stars per pc 2 , we identified six new clusters of YSO 
in the complex and a group located south from S256 
that may be part of the cluster in S256 (see Fig. 12). 
The first new cluster (G192.75-0.00), is located south of 
S258. The second one, called G192. 75-0.08, is located 
south-west from G192.75-0.00. The third (G192.63-0.00) 
is located east from S255-2, is more dispersed than the 
other two and includes the ionizing source of the Hn re- 
gion S255. The fourth cluster called G192.65-0.25, is lo- 
cated south-west from S256. The fifth and sixth clusters 
(G192.65-0.08 and G192.55-0.01) are located between the 
Hii regions S255 and S257, south from S255-2 and north 
from S255N respectively. The number of members in 
each cluster, as well as the number of Class I and Class 
II sources per cluster are shown in Table 4. We also iden- 
tified the already known clusters S255-2, S256, S258, and 
G192. 54-00. 15, located inside Hii region S254 (Dutra & 
Bica 2001). Finally, we calculated that the ratio of young 
stars in clusters vs. the total number of YSO is 0.8. If 
we change the value of d c in ±20%, the ratio change to 
0.83 and 0.75 respectively. 

3.4. Cloud structure 
3.4.1. K-band extinction map 

The K-band extinction map of the complex is shown in 
Figure 14. The map was created using the mean value of 
the K-band extinction of the nearest N = 10 stars over 
a uniform 3 arc-seconds grid. The Ak values were esti- 
mated by dereddening the sources with no IR-excess in 
the H-K vs. J-H color-color diagram, using the reddening 
law by Indebetouw et al. (2005). Finally, we convolved 
the resulting map with a 30 arc-seconds Gaussian kernel. 
The total number of sources used to build the map was 
1937. 

Since the complex is located at 2.4 kpc, foreground 
field stars will be reddened by interstellar extinction. The 
histogram in Figure 14 shows the distribution of Ak in 
the complex. If there were no cloud in the field, the 
value of Ak would rise, and then fall when it reaches our 
completeness limit, having a single peak. The presence 
of the molecular cloud splits the distribution of Ak in 
two groups since the value of the extinction rises where 
the cloud is located, separating the foreground from the 
background stars. This separation is seen in the his- 
togram at Ak ~ 0.4. This value is in good agreement 
with values of interstellar extinction at 2.4 kpc (0.15 mag 
per kpc in K-band, Indebetouw et al. (2005)). Therefore, 
to eliminate the foreground contribution to the extinc- 
tion, we only used stars with Ak > 0.4 in our extinction 
map. Also, since we are interested in the extinction due 
to only the molecular cloud (and not due to the ISM), 
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we shifted A^- down by 0.4 magnitudes. 

The Ak values derived may be under-estimated due to 
the small number of background stars detected. The 90% 
completeness limit in the K-band for our dataset is 16 
magnitudes, corresponding to a K0 type star at 2.4 kpc 
from the Sun. For example, if we want to achieve an 
extinction of 3 magnitudes in the K-band, we have to 
detect a star of spectral type AO at 2.4 kpc, or B8 at 3.0 
kpc. In this mass range, the number of background stars 
available to measure high extinctions decreases rapidly. 

The contours shown in Figure 14 start at Ak of 0.4 
and increase by increments of 0.2. The extinction map 
traces what appear to be three dusty filaments in the 
complex; extending north-south and centered on S255- 
2; south of S254 and S256, extending east-west; and 
from S255-2 to the east, in the middle of the Hn re- 
gion S255, where G192.63-0.00 is located. There are 
also regions with high extinction towards S258, G192.54- 
0.15, G192. 75-0.00, and G192. 75-0.08, and south from 
G192.69-0.25. Although the extinction map resembles 
the YSO density map presented in Figure 12, we notice 
one difference. The peaks in the extinction and in the 
YSO density maps for S256 and S258 are separated by 
approximately 1 arc-minute (0.6 pc at a distance of 2.4 
kpc) , while for the other groups of stars the peaks seem to 
match in position (see Figure 15). This difference could 
be explained by the fact that the Hn regions are remov- 
ing the dust from the clusters (see comparison with 13 CO 
column density in next section). 

3.4.2. Molecular Gas 

We used the column density of 13 CO and the velocity 
maps of 12 CO to study the structure of the molecular 
cloud. We derived the 13 CO column density from: 

A^ 3 CO) = 2.42 x 10" ( y^S^t ) cm" 2 , (3) 

where 

T « = 7 '"'La A ' ^ 

In (1 + Ta+0 S2 j 
T = -ln\l- 5 J 9 ex - j , (5) 

\ e 5 - 29 / T e*_l 1 / 

Ta is the observed peak temperature of the 12 CO, 
and Av correspond to the full width at half maximum 
(FWHM) of the 13 CO line (Heyer et al. 1989). We ob- 
tained the FWHM by fitting a Gaussian profile to the 
velocity spectrum in each pixel 8 . 

The 13 CO column density map is shown in Figure 16. 
There is a long (20pc) horizontal filament with a verti- 
cal component between S255 and S257. This filament is 
connected on the east side with the region S255B to the 
south (Bica et al. 2003), and with another vertical fila- 
ment to the north. There is also more molecular material 
south from the S254-S258 complex. The expanded view 
in Figure 16 shows the column density in the same area 
as the K-band extinction map, Hn regions and clusters 
arc indicated in this figure. The column density peak 
is located between Hn regions S255 and S257, forming 
a vertical filament in the middle of the Hn regions. We 

8 Based on code by R. Franco-Hernandez 



detected molecular material in the locations of all the 
clusters in the complex. Compared with the map from 
Heyer et al. (1989), the column density values obtained 
are similar. In our case, we do detect molecular mate- 
rial in the location of cluster G192.75-0.00. The same for 
cluster G192.54-0.15, but with a lower column density. 

In general the 13 CO density map agrees well with the 
K-band extinction map, showing the presence of gas 
around the locations where star formation is occurring, 
as well as in places without clusters but that are part 
of the molecular cloud, like south of S256 and west of 
G192.75-0.08. The column density peaks agree reason- 
ably well with the extinction peaks for most of the clus- 
ters. In S255-2, the column density peak is located 1 
arc-minute south of the extinction peak, this difference 
may be due a sensitivity problem in our extinction map 
caused by the non detection of background stars in places 
of high extinction. 

In both the column density map as well as in the ex- 
tinction map, there is no material inside the Hn region 
S254 other than in the location of G192.54-0.15, sug- 
gesting that the Hn region already removed most of the 
molecular gas and dust. If this is the case, the molecular 
material in G192.54-0.15 must likely show some evidence 
of interaction with the ionizing star (like a cometary 
structure) due to their proximity. Such an interaction is 
not seen in the IRAC mosaics, suggesting that the cluster 
may not be related with the complex. The 12 CO velocity 
maps in Figure 17 show that the main body of the molec- 
ular cloud has a Vlsr between 5 and 9 km s _1 , while 
the 12 CO emission in G192.54-0.15 has a Vlsr of 24- 
25 km s _1 . Since it seems unlikely that the 12 CO emis- 
sion at that position is not related with G192. 54-00. 15, 
the difference in the velocity with respect to the main 
body of the molecular cloud suggests that the cluster is 
located in the background of the complex (see § 4.1.4). 

Finally, we used the 13 CO column density to estimate 
the gas mass for clusters in the complex. In order to do 
this, we enclosed the cluster members by an ellipsoid and 
then we calculated the mass inside the ellipsoid assum- 
ing an abundance of 12 CO to 13 CO of 50 (Israel et al. 
2003), and of 12 CO to H 2 of 8.5xl0~ 5 (Frerking et al. 
1982). Gas masses, and other properties of the clusters 
are shown in Table 4 and are discussed in § 4. 

4. DISCUSSION 

4.1. Clusters in the complex 

We detected in total 510 sources with IR excess in the 
complex S254-S258, including 87 Class I, and 165 Class 
II sources. Of those, 398 are located in clusters (as de- 
scribed in § 3.3.2). Here, we discuss the global properties 
of the clusters. 

4.1.1. S255-2 & S255N 

These are located between Hn regions S255 and S257. 
We detected 23 Class I, 14 Class II, and 140 sources with 
IR-excess. We are likely to be missing some YSO due to 
the bright PAH emission (see § 3.2.2). The core of S255- 
2 is saturated in the IRAC bands. The peak surface 
density of sources with IR-excess is 156 stars pc~ 2 , and 
the K-band extinction peak is 1.4 magnitudes, located 
1 arc-minute south of S255-2. The peaks are coincident 
in location within 1 arc-minute (0.6 pc at a distance of 
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2.4 kpc). The 12 CO spectrum peaks at 7.5 km s 1 for 
S255-2 and at 8.4 km s" 1 for S255N. 

The detection of ultra-compact Hn regions (UCHu), 
outflows, and the high ratio of Class I to Class II suggest 
that S255-2 and S255N arc clusters in an early stage of 
formation (see also Cyganowski et al. (2007)). 

4.1.2. S256 & S256 south 

The YSO in S256 are more concentrated compared 
to S256-south. We detect 14 Class I, 44 Class II, and 
115 sources with IR-excess. The surface density peak is 
256 stars pc~ 2 and the peak in Ak is 1.3 magnitudes, 
located 1 arc-minute south-west from the surface den- 
sity peak. The 12 CO spectrum shows a double peak in 
S256, at 6.1 and 8.5 km s _1 , and a small component 
at 3.1 km s _1 . For S256-south there is a single peak 
at 7.4 km s _1 , and a small peak at 3.3 km s" 1 . The 
spectral type of the central source in the Hn region S256 
derived from our data is B0.9V. The source is classified 
as a Class II. 

Figure 18 shows 12 CO position-velocity maps of the 
area around S256. The peak at 8.5 km s _1 corresponds 
to a blob of molecular material that seems to have been 
ejected from the cloud. The mass of the blob is 30 M . 
If we assume a velocity of 2 km s^ 1 with respect to the 
main body of the cloud, the estimated kinetic energy is 
lxlO 45 ergs (the energy radiated by an B0.9V star in 
lxlO 5 yr is about lxlO 46 ergs). 

The presence of an Hii region instead of an UCHu re- 
gion, and also the ratio of Class I to Class II sources 
suggest that S256 and S256-south are more evolved than 
S255-2. 

4.1.3. S258 

We detect 7 Class I, 9 Class II, and 49 sources with 
IR-excess. The peak surface density is 298 stars pc~ 2 
and the peak Ak is 1.0 magnitudes. The 12 CO spec- 
trum has a strong peak at 8.0 km s _1 , and a small peak 
3.8 km s _1 . The spectral type of the central source in the 
Hii region derived from our data is B1.5, corresponding 
to a Class II. 

Similar to cluster S256, S258 seems more evolved than 
S255-2. The location of the central source in the Hn re- 
gion indicates that it is likely the ionizing source. The 
13 CO column density map shows that S258 is part of the 
same filament where S255-2 and S256 are located. 

4.1.4. G192.54-0.15 

This cluster is located inside the Hii region S254. We 
detect 2 Class I, Class II, and 9 sources with IR- 
excess. The peak in density is 166 stars pc~ 2 and the 
peak Ak is 0.7 magnitudes. The 12 CO spectrum peaks 
at 24.5 km s _1 . Since the cluster is located near the 
center of the Hn region S254, where most of the molecu- 
lar material has been evacuated, we are able to detect it 
through the molecular cloud. 

The kinematic distance for G192.54-0.15, using the ro- 
tational model from Brand & Blitz (1993) is about 9 
kpc (the uncertainty in the estimation of the kinematic 
distance toward the anti-center of the Galaxy may be as 
much as a factor of 2, Reid (2008)). At that distance, the 
cluster will be located in the outer parts of the Galaxy, 
having an IRAC luminosity 4 times less than the trapez- 
ium in Orion. There is evidence of molecular material in 



the same direction with similar Vlsr on a larger scale 
(Carpenter et al. 1995). G192.54-0.15 may be part of 
star formation occurring in those molecular clouds, and 
it could be an important object for further study of star 
formation in the outer Galaxy. 

4.1.5. G 192.75-0.00 

This cluster is located to the south of S258. We detect 
Class I and 7 Class II sources. The peak in surface 
density is 76 stars pc~ 2 and the peak in Ak is 0.7 mag- 
nitudes. The 12 CO spectrum peaks at 8.1 km s -1 . 

The absence of Class I objects suggests that this cluster 
may be older than S256 and S258. There is no signature 
of an Hii region, implying that this cluster contains low 
mass stars. 

4.1.6. G 192.75-0.08 

Located south-east of G192. 75-0.00. We detect 25 
sources with IR-excess, 9 Class I and 15 Class II. The 
peak in surface density is 238 stars pc -2 and the peak 
in Ak is 1.1 magnitudes. The 12 CO spectrum peaks at 
7.2 km s _1 . 

As for G192. 75-0.00, in G192. 75-0.08 there is no sig- 
nature of massive star formation. But the presence of 
Class I sources suggests that it is younger than G192.75- 
0.00. The value of Ak, as well as the surface density of 
gas, is higher than in G192. 75-0.00. 

4.1.7. G 192. 63- 0.00 

This cluster is located inside Hn region S255. We de- 
tect 7 Class I, 5 Class II and 22 sources with IR-excess. 
The peak in surface density is 39 stars pc~ 2 and the peak 
in Ak is 1.0 magnitudes. The 12 CO spectrum peaks at 
6.9 km s _1 . The ionizing star of S255 is located inside 
this cluster. The spectral type of the ionizing star derived 
from our data is B0.0, and is classified as a Class II. 

The IRAC classification of the ionizing star in S255 
(Class II) suggests that it is in an earlier stage than the 
ionizing star in S257 (Class III). In this context, we may 
assume that other YSO formed together with the ionizing 
star of S255 are also in an earlier stage than the YSO 
formed with the ionizing star of S257. In this case, we 
would expect to detect more sources with IR-excess in 
S255 than in S257. In fact, we detect almost 3 times 
more IR-excess sources inside Hn region S255 than in 
S257. We hypothesize that G192. 63-0.00 is the natal 
cluster of the ionizing star in S255. 

4.1.8. G192.69-0.25 

This cluster is located south-west from S256. We de- 
tect 7 sources with IR-excess; 3 Class I and 4 Class II. 
The peak in surface density is 27 stars pc -2 and the peak 
in Ak is 0.8 magnitudes. The 12 CO spectrum peaks at 
7.0 km s" 1 . 

Since this cluster is located at the edge of our mosaic, 
we may be missing some members. The 13 CO column 
density map shows more molecular material in the loca- 
tion of the cluster south of our mosaic. 

4.1.9. G192.65-0.08 & G192.55-0.01 

These two clusters are located between Hn regions 
S255 and S257. G192.65-0.08 is located south from S255- 
2. We detect 13 sources with IR-excess; 2 Class I and 6 
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Class II. The peak in surface density is 37 stars pc~ 2 and 
the peak in Ak is 0.9 magnitudes. The 12 CO spectrum 
peaks at 7.0 km s _1 . G192.55-0.01 is located north from 
S255N. We detect 9 sources with IR-excess; 2 Class I and 
5 Class II. The peak in surface density is 30 stars pc~ 2 
and the peak in Ak is 0.7 magnitudes. The 12 CO spec- 
trum peaks at 6.3 km s _1 . 

4.2. Star formation efficiency 

In order to estimate the star formation efficiency of 
clusters: 



Stromgren sphere given by Stahler & Palla (2004), 



M, 



e = 



+ M, 



(6) 



gas 



we need to estimate the mass of gas and stars in the 
clusters. The mass of gas was calculated using the 13 CO 
column density (see § 3.4.2). To estimate the mass of 
stars in clusters we count the number of stars and assume 
each one has a mass of 0.5 M Q . To count the number of 
stars we must correct for two sources of incompleteness: 
1) cluster members without IR-excess are not identified 
by our method, and 2) we do not detect cluster members 
fainter than 16 magnitudes in the K-band. To account 
for cluster members without IR-excess we calculated the 
K-band surface density in a region of the map without 
IR-excess sources (control field), then we calculated the 
number of cluster members by subtracting the mean sur- 
face density (normalized by the area) of the control field 
from the mean surface density of the clusters. To es- 
timate the number of cluster members fainter than 16 
magnitudes in K-band, we created a 1 Myr old synthetic 
cluster having an initial mass function (IMF) given by 
the Orion nebula cluster. We compared the observed 
K-band luminosity function (KLF) of a region including 
S255-2, S255N, S256 and S258 with the resulting KLF 
of synthetic clusters with ages of 0.5, 1, 3, and 5 Myr 
(Muench et al. 2000). We found that the 1 Myr old clus- 
ter best represents the young population in those regions. 
Then, using the synthetic 1 Myr old cluster, we estimated 
the percentage of cluster members brighter than our com- 
pleteness limit (16 magnitudes in K-band), at a distance 
of 2.4 kpc and assuming an extinction of 0.4 magnitudes 
in the K-band. We found that this percentage is around 
50%, therefore the corrected number of cluster members 
is two times the number of observed members. We then 
estimated the mass of stars in clusters by multiplying 
the corrected number of cluster members by 0.5, which 
corresponds to the average mass of stars. 

Star formation efficiencies for clusters in the complex 
are shown in the last column of Table 4. We estimate 
an uncertainty of a factor of 2 in the values of the ef- 
ficiency. The values are less than 20% and comparable 
with those in other regions of star formation (Elmegreen 
et al. 2000) for most of the clusters with the exception 
of G192.54-0.15, G192. 75-0.00, and G192.55-0.01. The 
cluster G192.54-0.15 has an efficiency of 54% at 2.4 kpc 
from the Sun. The efficiency of this cluster at 9 kpc from 
the Sun is 1%. 

4.3. Ages of the Hn regions and relative ages of clusters 

We considered the expansion time of the Hn re- 
gions using the relation for the second expansion of the 



t = 



4rs 
7c 



7/4 



where the Stromgren radius is: 



( 37V 
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r is the size of the Hn region, c is the velocity of the 
sound, and a is the recombination coefficient (Pcquignot 
et al. 1991). The sizes of the Hn regions were estimated 
using Ha imaging from Mizuno (1982). We assumed a 
Hydrogen density of njj = 3 x 10 4 cm~ 3 (Wright et al. 
1981), a gas temperature of 10 4 K and a sound speed of 
11 km s -1 . We use the spectral type given by our data 
(listed in Table 3) to estimate the amount of ionizing pho- 
tons Af. The estimated ages for the Hn regions in ascend- 
ing order are: S258 (lxl0 5 yr), S256 (2xl0 5 yr), S255 
(1.5xl0 6 yr), S257 (1.6xl0 6 yr), and S254 (5.0xl0 6 yr)- 

We used the ratio of Class I to Class II sources (I/II) 
to estimate the relative ages of clusters in the complex. 
If we assume that all the stars in a cluster were formed at 
approximately the same time, and since Class I objects 
arc in an earlier stage of formation than Class II objects, 
we expect a high I/II ratio in young clusters. We also 
used the ratio of Class I to IR-excess sources (I/N IR ) to 
probe clusters with small number of members. The ratios 
I/II and I/N IR per cluster are shown in Table 4. Even 
though the relative ages for the clusters in the complex 
given by both ratios is different, they give the same or- 
der from youngest to oldest in clusters S255-2, S258 and 
S256. 

4.4. Star formation in S254-S258: the overall picture 

We estimated the ages of the ionizing stars in the com- 
plex using the expansion time of their Stromgren sphere. 
From young to old the Hn regions are ordered as follows: 
S258, S256, S255, S257, and S254. Since the spectral 
type of the ionizing star in S255 and S257 is similar (the 
same for S256 and S258, see Table 3), we can assume that 
both stars are evolving with similar time-scales. In this 
context, the age estimated using the Stromgren sphere is 
consistent with their IRAC+NIR classification (Class II 
for S255 and Class III for S257). We found the same for 
the ionizing stars in S256 (Class II) and S258 (Class I). 
Also, the relative ages estimated for clusters S256 and 
S258 using the ratios I/II and I/N IR (S258 relatively 
younger than S256) is consistent with the ages estimated 
from the Hn regions. 

Both the column density map and the extinction map 
show that most of the gas and dust inside the Hn region 
S254 has been removed. For Hn regions S255 and S257, 
there is still some molecular material inside, while S256 
and S258 are still embedded in the molecular cloud. This 
sequence is consistent with the evolution of the complex 
suggested by the ages of the Hn regions. The column 
density map and the extinction map show that a filament 
of material is still present inside the Hn region S255, 
while in S257 there is less molecular material. This is 
expected if S257 is older than S255 (since the spectral 
type of the ionizing stars is similar, the evolution of the 
Hn regions should also be similar). This is consistent 
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with the difference in the 1RAC+NIR classification of 
the ionizing sources (S255 in an earlier stage than S257), 
and with the detection of 3 times more IR-excess sources 
inside S255 than in S257. Those sources may be part of 
the cluster where the ionizing star of S255 was born. 

The K-band extinction map shows that the extinction 
in clusters S256 and S258 peaks in a different location 
than the surface density of YSO. This is similar to Orion, 
where part of the star forming region has been exposed 
by the action of the massive stars. In the case of S255-2, 
the Hn region is not formed yet, and for G192.75-0.00, 
G192. 75-0.08, and G192.69-0.25, there is not evidence of 
massive stars in those clusters. 

Using the 12 CO velocity maps in Figure 17 and the 
IRAC mosaics, we can study the spatial distribution of 
the Hn regions in the complex. The molecular material 
accumulated between Hn regions S255 and S257 suggests 
that those Hn regions are located at the same distance 
from the Sun. For S254, the absence of accumulated ma- 
terial between S257 and S254 or any evidence of interac- 
tion with S257 or evidence that star formation happened 
between S254 and S257 (like in between S255 and S257), 
suggests that it is located outside the plane of S255 and 

5257. Finally, the relative Vlsr peaks between G192.63- 
0.00, S255-2', and S255N (6.9, 7.5, and 8.4 km s" 1 respec- 
tively) suggests that star formation is propagating to- 
ward the back of the Hn regions S255 and S257, (the rel- 
ative kinematic distance between those clusters is more 
than 150 pc). 

5. CONCLUSIONS 

We used Spitzer-IRAC and Flamingos-NIR data to 
identify YSO in the massive star forming complex S254- 

5258. We used the IRAC slope of the SED a IKAC to 
classify YSO into Class I and Class II sources, and the 
K-4.5 /an v/s H-K color-color diagram to identify sources 
with IR-excess. We detected in total 87 Class I, and 165 
Class II sources. Also, we identified in total 510 sources 
with IR-excess. 

We used the surface density map of sources with IR- 
excess, and NN analysis to study the spatial distribution 
of YSO. We discovered six unknown clusters of YSO in 
the complex: G192.63-0.00, G192.75-0.00, G192. 75-0.08, 
G192.69-0.25, and clusters G192.65-0.08 and G192.55- 
0.01 that are located between Hn regions S257 and S255. 
We also discovered a southern component of the cluster 
S256. We found that 80% of the IR-excess sources are 
located in clusters, while the rest is distributed in a more 
isolated component. 
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We estimated the ages of the ionizing stars in the 
complex using the expansion time of their Stromgrcn 
sphere. The result is consistent with the IRAC+NIR 
classification of the ionizing stars, and with the rela- 
tive ages estimated using the I/II and I/N IR ratios. The 
ages estimated for the Hn region are: S258 (lxl0 5 yr), 
S256 (2xl0 5 yr), S255 (1.5xl0 6 yr), S257 (1.6xl0 6 yr), 
and S254 (5.0x 10 6 yr)- 

We used the 13 CO column density and the Ak extinc- 
tion maps to study the structure of the molecular cloud. 
The main body of the molecular cloud has Vlsr between 
5-9 km s _1 . The Hn regions S255 and S257 are located 
in the same plane while S254 seems to be located behind 
S255 and S257. 

The Vlsr of cluster G192.54-0.15 suggests that it 
is located in the background of the complex. The 
double peak in the 12 CO spectrum of S256 and 12 CO 
position-velocity maps of S256 suggests that a blob of 
molecular material is being pushed to the back of the 
molecular cloud by the Hn region. 
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TABLE 1 

Summary of detections' 1 and areas observed per filter 





J 


H 


K 


H b 


K b 


JHK C 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


All IRAC 


HK[4.5] C 


Detections 


5941 


5966 


7430 


152 


145 


4225 


6131 


6391 


1472 


853 


462 


3838 


Limiting magnitude d 


16.0 


16.4 


16.0 


13.6 


13.8 




14.0 


13.5 


11.0 


10.0 






Area [arc-min 2 ] 


530 


530 


530 


1 


1 


530 


500 


570 


500 


570 


430 


500 



a Detections with error < 0.2 magnitudes. 
b NIRC data. 

c Including Flamingos and NIRC data. 
d Magnitudc for 90% completeness. 
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TABLE 2 
Photometry 11 



RA (2000) Dec (2000) 



K 



[3.6] 



[4.5] 



[5.8] 



err Class 



93 


36916 


17.92738 


14 


43 





01 


12 


91 





01 


12 


02 


0.01 


11 


57 


0.04 


11 


17 


0.04 


10 


04 


0.13 


8 


26 





13 


93 


28009 


17.98375 


17 


91 


o 


08 


15 


78 





03 


14 


16 


0.01 


12 


09 


0.02 


10 


94 


0.02 


9 


76 


0.02 


8 


59 





03 


93 


10282 


17.98915 


15 


86 


o 


03 


14 


71 


o 


02 


13 


65 


0.02 


12 


39 


0.07 


11 


79 


0.05 


10 


21 


0.15 


8 


68 


o 


19 


93 


14838 


17.94358 


15 


72 


o 


01 


13 


90 


o 


01 


12 


77 


0.01 


11 


12 


0.02 


10 


35 


0.01 


9 


85 


0.04 


9 


01 


o 


09 


93 


15528 


17.94797 


15 


46 


o 


02 


14 


03 


o 


01 


13 


08 


0.01 


11 


71 


0.02 


11 


07 


0.02 


10 


74 


0.05 


9 


61 


o 


08 


93 


15188 


17.94479 


16 


17 


o 


04 


14 


22 


o 


01 


13 


04 


0.01 


11 


85 


0.02 


11 


29 


0.02 


10 


81 


0.05 


9 


87 


o 


11 


93 


24628 


18.04524 


15 


37 


o 


01 


14 


13 


o 


01 


13 


32 


0.01 


12 


19 


0.02 


11 


74 


0.02 


11 


16 


0.03 


9 


91 


o 


03 


93 


15045 


17.94815 


14 


27 


o 


01 


13 


18 


o 


01 


12 


53 


0.01 


11 


51 


0.02 


11 


09 


0.02 


10 


94 


0.04 


9 


96 


o 


08 


93 


21329 


18.01470 


18 


37 


o 


11 


15 


82 


o 


02 


14 


21 


0.01 


12 


35 


0.02 


11 


90 


0.02 


11 


06 


0.05 


10 


02 


o 


09 


93 


23944 


18.02592 


17 


47 


o 


06 


16 


02 


o 


03 


14 


76 


0.02 


13 


61 


0.04 


12 


34 


0.03 


11 


67 


0.04 


10 


12 


o 


05 


93 


35432 


17.88051 


16 


42 


o 


03 


14 


79 


o 


01 


13 


64 


0.01 


12 


02 


0.02 


11 


22 


0.02 


10 


99 


0.02 


10 


12 


o 


02 


93 


16369 


17.94083 


16 


80 


o 


03 


14 


95 


o 


01 


14 


13 


0.01 


12 


84 


0.04 


12 


35 


0.03 


11 


59 


0.11 


10 


12 


o 


15 


93 


25438 


17.96418 


17 


38 


o 


06 


15 


46 


o 


02 


14 


14 


0.01 


12 


38 


0.03 


11 


71 


0.02 


11 


21 


0.07 


10 


25 


o 


13 


93 


36552 


17.94602 


16 


45 


o 


02 


14 


93 


o 


01 


13 


94 


0.01 


12 


67 


0.03 


12 


14 


0.02 


11 


51 


0.04 


10 


26 


o 


06 


93 


07682 


17.98306 


16 


08 


o 


01 


15 


18 


o 


01 


14 


70 


0.01 


13 


74 


0.04 


13 


01 


0.04 


12 


18 


0.04 


10 


28 


o 


03 


93 


11366 


17.98991 


17 


55 


o 


04 


16 


29 


o 


03 


15 


21 


0.01 


13 


81 


0.04 


13 


36 


0.04 


12 


05 


0.05 


10 


46 


o 


05 


93 


38181 


17.93109 


17 


09 


o 


04 


15 


70 


o 


03 


14 


75 


0.01 


13 


80 


0.06 


13 


47 


0.05 


12 


36 


0.12 


10 


54 


o 


15 


93 


16968 


18.02927 


16 


04 


o 


01 


15 


04 


o 


01 


14 


56 


0.01 


13 


80 


0.05 


13 


63 


0.05 


12 


15 


0.09 


10 


56 


o 


10 


93 


25941 


18 00423 


15 


62 


o 


01 


14 


42 


o 


01 


13 


68 


0.01 


12 


95 


0.03 


12 


59 


0.03 


11 


77 


0.05 


10 


60 


o 


07 


93 


22853 


18.04282 


16 


54 


o 


02 


15 


42 


o 


02 


14 


60 


0.01 


13 


35 


0.03 


12 


70 


0.03 


11 


83 


0.04 


10 


63 


o 


06 


93 


15380 


17.94118 


17 


07 


o 


04 


15 


55 


o 


02 


14 


73 


0.02 


13 


46 


0.05 


12 


72 


0.03 


11 


82 


0.09 


10 


68 


o 


15 


93 


26565 


17.95645 


16 


58 


o 


03 


15 


08 


o 


01 


14 


42 


0.01 


13 


60 


0.04 


13 


43 


0.05 


12 


29 


0.08 


10 


74 


o 


08 


93 


27636 


17.98397 


16 


16 


o 


02 


14 


45 


o 


01 


13 


43 


0.01 


12 


15 


0.02 


11 


62 


0.02 


11 


35 


0.05 


10 


76 


o 


12 


93 


12441 


17.91705 


18 


71 


o 


10 


15 


58 


o 


01 


13 


93 


0.01 


12 


36 


0.02 


11 


99 


0.02 


11 


12 


0.02 


10 


79 


o 


03 


93 


22486 


18.01707 


15 


00 


o 


01 


14 


48 


o 


01 


13 


91 


0.01 


12 


67 


0.03 


12 


01 


0.02 


11 


11 


0.08 


10 


80 


o 


14 


93 


24936 


17.99069 


16 


45 


o 


02 


15 


15 


o 


01 


14 


57 


0.01 


13 


46 


0.04 


13 


03 


0.04 


11 


83 


0.14 


10 


80 


o 


17 


93 


07597 


17.87753 


16 


77 


o 


03 


14 


66 





01 


13 


62 


0.01 


12 


63 


0.02 


12 


12 


0.02 


11 


36 


0.02 


10 


83 





03 


93 


27975 


17.80454 


16 


30 


o 


02 


14 


91 


o 


01 


13 


88 


0.01 


12 


73 


0.03 


12 


32 


0.03 


11 


46 


0.02 


10 


84 


o 


03 


93 


23428 


18.04390 


16 


60 


o 


02 


15 


59 


o 


02 


14 


82 


0.01 


13 


86 


0.04 


13 


26 


0.04 


12 


74 


0.07 


11 


00 


o 


07 


93 


11582 


17.91519 


18 


00 


o 


05 


16 


11 


o 


02 


15 


06 


0.01 


13 


40 


0.03 


13 


24 


0.04 


11 


82 


0.03 


11 


08 


o 


04 


93 


13181 


17.93129 


15 


60 





01 


14 


46 





01 


14 


03 


0.01 


13 


34 


0.03 


12 


98 


0.04 


12 


39 


0.08 


11 


10 





10 


93 


12912 


17.91860 


18 


47 


o 


09 


16 


15 


o 


02 


14 


73 


0.01 


13 


62 


0.04 


12 


91 


0.03 


11 


80 


0.03 


11 


15 


o 


03 


93 


29075 


18.04884 


15 


69 


o 


03 


14 


95 


o 


02 


14 


20 


0.01 


13 


37 


0.03 


12 


84 


0.03 


12 


29 


0.04 


11 


16 


o 


05 


93 


20743 


18.07821 


15 


33 





02 


14 


17 





01 


13 


41 


0.01 


12 


60 


0.02 


12 


24 


0.02 


11 


79 


0.03 


11 


18 





04 


93 


34914 


17.88210 


16 


22 





02 


14 


61 





01 


13 


71 


0.01 


12 


96 


0.03 


13 


26 


0.04 


11 


60 


0.02 


11 


23 





03 


93 


21317 


18.02341 


16 


64 





06 


15 


59 





04 


14 


65 


0.02 


13 


47 


0.04 


13 


04 


0.04 


12 


36 


0.09 


11 


27 





16 


93 


09181 


17.84300 


15 


71 





01 


14 


52 





01 


13 


90 


0.01 


13 


18 


0.03 


12 


81 


0.03 


12 


32 


0.04 


11 


30 





03 


93 


03669 


17.81515 


16 


58 





03 


15 


38 





02 


14 


48 


0.02 


13 


47 


0.04 


12 


54 


0.03 


12 


46 


0.04 


11 


31 





03 


93 


28318 


18.05870 


16 


98 





05 


16 


44 





05 


15 


20 


0.02 


13 


85 


0.04 


13 


22 


0.04 


12 


61 


0.05 


11 


37 





04 


93 


18279 


17.94522 


16 


60 





02 


15 


00 





01 


14 


12 


0.01 


12 


94 


0.03 


12 


60 


0.03 


12 


09 


0.05 


11 


39 





10 


93 


09937 


17.95723 


15 


34 





01 


14 


33 





01 


13 


79 


0.01 


13 


02 


0.03 


12 


32 


0.03 


12 


14 


0.06 


11 


40 





13 


93 


26479 


18.03973 


17 


04 





03 


16 


04 





03 


15 


18 


0.02 


14 


04 


0.05 


13 


53 


0.05 


12 


68 


0.10 


11 


40 





17 


93 


05625 


17.98123 


15 


76 





01 


14 


86 





01 


14 


21 


0.01 


13 


44 


0.03 


12 


89 


0.03 


12 


44 


0.04 


11 


45 





04 


93 


08435 


17.81195 


16 


03 





01 


14 


71 





01 


14 


12 


0.01 


13 


48 


0.04 


12 


86 


0.03 


12 


30 


0.04 


11 


46 





03 


93 


27053 


17.85008 


19 


31 





22 


16 


69 





05 


15 


18 


0.02 


13 


87 


0.04 


13 


47 


0.04 


12 


16 


0.03 


11 


47 





04 


93 


12422 


17.93572 


15 


90 





01 


14 


74 





01 


14 


23 


0.01 


13 


21 


0.03 


13 


03 


0.04 


12 


28 


0.04 


11 


49 





07 



a The entire table is available electronically, in the on-line version of this paper. 
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TABLE 3 
Parameters of ionizing stars 



Name 


Name 


RA (2000) 


Dec (2000) 


Spectral Typc a 


Spectral Type 


HII region age 


Reference 




Catalog 


hh mm ss 


dd mm ss 


Literature 


Our data 


[Yr] 






S254 


HD253247 


06 12 22.1 


18 00 57.7 


O9.0V 


09.6V 


5.1xl0 6 


Moffat et al. 


(1979) 


S255 


HD253237 


06 13 04.2 


17 58 41.1 


B0.0V 


B0.0V 


1.5xl0 6 


Moffat et al. 


(1979) 


S256 




06 12 36.5 


17 56 54.3 


B2.5V 


B0.9V 


2xl0 5 


Russeil et al. 


(2007) 


S257 


ALS19 


06 12 44.2 


17 59 13.8 


B0.5V 




1.6xl0 6 


Moffat et al. 


(1979) 


S258 




06 13 28.3 


17 55 33.3 




B1.5V 


1 x 10 5 







a Spectral types from literature correspond to the reference listed in the last column of the table. 
b Agcs were calculated using the spectral type given by our data. For S257 we used the spectral type from Moffat 
et al. (1979). 
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TABLE 4 
Parameters of clusters 21 



Name RA (2000) 
hh mm ss 


Dec (2000) 
dd mm ss 




N c 


I 


II 


I/II 


I/N IR 


0"s 

[pc 


d 

- 2 ] 


cr 3m 
[pc" 


e 
-2] 


[cm" 2 ] 


rr 8 

' ' gmax 

[cm" 2 ] 


[M ] 


e 


S255-2 & S255N 


6 


12 


54 


17 


59 


10 


140 


488 


23 


14 


1.64 h 


0.16 




18 




156 


2.9xl0 16 


7.0xl0 16 


1890 


0.11 


S256 


6 


12 


36 


17 


56 


48 


115 


276 


14 


44 


0.32 


0.12 




16 




256 


1.6xl0 16 


3.6xl0 16 


2005 


0.06 


S258 


6 


13 


29 


17 


55 


44 


49 


186 


7 


9 


0.78 


0.14 




17 




298 


1.4xl0 16 


2.3xl0 16 


815 


0.10 


G192.54-0.15 


6 


12 


24 


17 


59 


27 


9 


88 


2 







0.22 




17 




166 


2.0 xlO 15 


4.1xl0 15 


40 


0.54 


G192.75-0.00 


6 


13 


24 


17 


52 


51 


7 


60 





7 


0.00 


0.00 




12 




76 


8.0 xlO 15 


9.2xl0 15 


100 


0.23 


G192. 75-0.08 


6 


13 


05 


17 


50 


30 


25 


50 


9 


15 


0.60 


0.36 




11 




238 


1.2 xlO 16 


1.9xl0 16 


490 


0.05 


G192.63-0.00 


6 


13 


08 


17 


58 


42 


22 


44 


7 


5 


1.40 h 


0.32 




10 




39 


1.4xl0 16 


2.0xl0 16 


560 


0.04 


G192.69-0.25 


6 


12 


19 


17 


48 


28 


7 


14 


3 


4 


0.75 


0.43 




9 




27 


l.OxlO 16 


1.4xl0 16 


125 


0.05 


G192.65-0.08 


6 


12 


53 


17 


55 


30 


13 


26 


2 


6 


0.33 


0.15 




9 




37 


l.SxlO 16 


3.0xl0 16 


340 


0.04 


G192. 55-0.01 


6 


12 


56 


18 


02 


44 


9 


18 


2 


5 


0.40 


0.22 




9 




30 


2.3 xlO 15 


4.8xl0 15 


20 


0.33 



a Wc assumed a distance of 2.4 kpc from the Sun in the estimation of physical parameters. 
^Number of stars with IR-excess, including Class I and Class II. 
c Corrected number of cluster members, see § 4.2. 
^Mcan surface density of stars. 
e Peak surface density of stars. 

Mean 13 CO surface density. 
s Peak 13 CO surface density. 

h Likcly to have missed some Class II (see end of § 3.2.2). 
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4Q-Q 30 20.0 1 0.0 8:13:00.0 50.O 4Q.0 3Q.0 ZD.O 12:10.0 



Fig. 1. — IRAC 3-color image of the S254-S258 complex (blue: 3.6/^m, green: 4.5/xm, red: 8.0/xm). White circles enclose the HII regions 
(Mizuno 1982). Clusters S255-2 and S255N are also indicated. Black rectangle shows the area covered by NIRC data. Known ionizing 
stars are indicated with arrows. 
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Fig. 2. — Apparent magnitude histograms of detections for all IRAC bands. Dashed area shows detections with error < 0.2 magnitudes. 
Since 3.6 and 4.5/^m bands are more sensitive to the photospheres of stars than 5.8 and 8.0/^m bands, the number of sources detected at 
5.8 and 8.0/im is less than at 3.6 and 4.5/im. Dashed lines indicate 90% completeness. 
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Fig. 4. — Spectral energy distribution (SED) of a Class I (left) and Class II (right) source in the S254-S258 complex. The observed fluxes 
in the JHK and all IRAC bands are shown by diamonds. Dereddencd fluxes are shown by crosses. Dashed line shows the fitted Oirac- 
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Fig. 5. — Distribution of spectral indeces of sources detected in all IRAC bands. Dashed lines separate Class I and Class II sources. 
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Fig. 6. — IRAC color-color diagram of observed magnitudes for all sources detected in all four IRAC bands. Sources detected only 
in the IRAC bands (shown in black) were classified using their observed Oirac value as Class I (o/bac > 0, shown as circles), Class II 
(—2 < Chirac < 0, shown as squares), or Class III (etmAC < ~2, shown as points). Sources that were also detected in JHK (shown in 
red), were classified using their dereddened Qirac value. Large dashed rectangle indicates the position of Class II sources from Allen et al. 
2004. The reddening vector is from Flaherty et al. 2007. 
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Fig. 7. — Upper-left: Distribution of sources detected inside (gray), and outside (white) areas with bright PAH emission (see areas 
marked in Figure 3c) at 5.8 fim. Upper-right: Distribution of Class I (gray), and Class II sources at 5.8 fim. Lower: The same but at 
8.0 fim. We don't detect faint sources inside the area of bright PAH emission. 
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Fig. 8. — Photometric uncertainty of sources detected at 5.8 /im {upper) and 8.0 (im (lower). Sources located inside areas with bright 
PAH emission are marked as crosses, sources located outside those areas are marked as points. Sources located outside areas with bright 
PAH emission extend to fainter magnitudes. 
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Fig. 9. — Flamingos H, K and IRAC 4.5 fim band color-color diagram for S254-S258. The black line near (0,0) shows the main sequence 
locus of a late type M dwarf star (Patten et al. 2006). The long arrow corresponds to the reddening vector (from Flaherty ct al. 2007). 
Stars with infrared excess, shown as diamonds, lay between the reddening vector and the dashed horizontal line. 
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Fig. 10. — Upper-left: Color-magnitude diagram for sources in S254-S258 detected in all IRAC bands. Circles show Class I sources, 
squares show Class II sources and points correspond to Class III sources. Dashed lines delimit background contaminated area. Upper- 
right: Distribution of galaxies from the IRAC Shallow Survey for the same color-magnitude diagram. Contours begin at 1.0 and increase 
by 1.0. Lower-left: Color-magnitude diagram for sources with H, K, and IRAC 4.5 (im detections. Sources with IR-excess are shown 
in diamonds. Dashed lines delimit background contaminated area. Lower-right: Distribution of Shallow Survey galaxies for the color- 
magnitude diagram. Contours begin at 1.0 and increase by 1.0. 
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Fig. 11. — Upper: Sources identified as Class I are shown marked with circles on the IRAC 4.5 (im mosaic of S254-S258. Lower: Sources 
identified as Class II (squares) on S254-S258. Large black circles correspond to the HII regions in the complex (from Mizuno (1982)) 
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Fig. 12. — Sources with IR-excess in the S254-S258 complex. Cluster members are shown in black and red points, the difference in color 
is just to separate the clusters. Isolated sources are shown in black plusses. Black contours represent the n = 5 nearest neighbor surface 
density of YSO. Contour levels are at 5 and 10 stars pc -2 . HII regions are shown in red color. New clusters are labeled in blue. White 
ovals in the center enclose clusters S255-2 and S255N. 
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Fig. 13. — Cumulative distribution of nearest-neighbor (NN) distances for sources with IR-excess in the S254-S258 complex. The 
low-density component starts at a distance of 0.01 degrees (or 0.4 pc at 2.4 kpc of distance). 
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Fig. 14. — Upper: K-band extinction map in the S254-S258 complex. Contours begin at 0.4 Ax and increase by 0.2. Labels are the same 
as in Figure 12. Lower: Distribution of for sources in the complex. Dashed line indicates the foreground extinction value adopted. 
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Fig. 15. — K-band extinction map in the S254-S258 complex, overlaid with contours of the surface density of sources with IR-excess. 
Contour levels are at 5, 20 and 80 stars pc -2 . Labels are the same as in Figure 12. Note the ~ 1 arc-minute offset between the extinction 
and stellar density peaks in S256 and S258. 
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Fig. 16. — Upper: Map of 13 CO column density derived from the ratio of 12 CO and 13 CO intensities. Contours are spaced logarithmically 
between 6 X 10 15 and 6 X 10 16 cm~ 2 . Lower: Close up of the boxed region having the same size as the K-band extinction map (shown in 
Figure 14). Labels are the same as in Figure 12. 
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Fig. 17.— 12 CO channel maps. The velocity (km s ) is indicated by the number in the top-right corner. Contour levels begin at 
3K, and increase by 3K. HII regions are indicated by black circles, the approximate size of the clusters is indicated by segmented ovals. 
Black line represents a scale bar of 2 pc. As shown also by Heyer et al. (1989), the bulk of the cloud is located at velocities between 5 and 
10 km s _1 . The only emission that matches the position of G192.54-0.15 is detected at 24.5 km s _1 (lower right). 
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Fig. 18. — 12 CO position and position-velocity maps in galactic coordinates of the S254-S258 complex. In a), we show the emission 
map at 7 km s~ 1 . The black ovals represent the size of the HII regions, b) and c) are position-velocity maps in latitude and longitude 
respectively across the dashed lines in a), c) shows a hole in the molecular cloud located at the same position as S256, and a blob of gas 
along the same line of sight at 9 km s _1 . 



